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Abstract 

Suppose an electric current I flows along a magnetic flux tube that has poloidal flux tp and radius 
a = a(z) where z is the axial position along the flux tube. This current creates a toroidal magnetic 
field B^. It is shown that, in such a case, non-linear, non-conservative J x B forces accelerate plasma 
axially from regions of small a to regions of large a and that this acceleration is proportional to 
dl 2 /dz. Thus, if a current-carrying flux tube is bulged at, say, z = and constricted at, say, 
z = ±/i, then plasma will be accelerated from z = ±h towards z = resulting in a situation similar 
to two water jets pointed at each other. The ingested plasma convects embedded, frozen-in toroidal 
magnetic flux from z = ±/i to z = 0. The counter-directed flows collide and stagnate at z = 
and in so doing (i) convert their translational kinetic energy into heat, (ii) increase the plasma 
density at z ~ 0, and (hi) increase the embedded toroidal flux density at z ~ 0. The increase in 
toroidal flux density at z ~ increases Bj, and hence increases the magnetic pinch force at z ~ 
and so causes a reduction of a(0). Thus, the flux tube develops an axially uniform cross-section, 
a decreased volume, an increased density, and an increased temperature. This model is proposed 
as a likely hypothesis for the long-standing mystery of why solar coronal loops are observed to be 
axially uniform, hot, and bright. It is furthermore argued that a small number of tail particles 
bouncing between the approaching counterstreaming plasma jets should be Fermi accelerated to 
extreme energies. Finally, analytic solution of the Grad-Shafranov equation predicts that a flux 
tube becomes axially uniform when the ingested plasma becomes hot and dense enough to have 
IHqukI '/ 'Bp ol = (/xoia(0)/V>) 2 /2; observed coronal loop parameters are in reasonable agreement 
with this relationship which is analogous to having (3 po i = 1 in a tokamak. 
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I. INTRODUCTION 



A long standing mystery in solar physics is why solar coronal loops typically have an ax- 
ially uniform cross-section [lj]; i.e., a filamentary shape. This issue has been made especially 
pressing by recent TRACE (Transition Region and Coronal Explorer ) spacecraft soft x-ray 
images which show a multitude of highly-defined axially uniform loops [2j; for example, see 
Fig. 1. Axial uniformity of flux tubes is also commonly observed in laboratory experiments, 
for example, in recent simulations of solar prominences jjj. 




FIG. 1: TRACE soft X-ray photo showing axial uniformity of coronal loops (image courtesy 
Lockheed Martin Space and Astrophysics Lab). As discussed in the text, the slight wrapping of 
flux tube #1 around flux tube #2 indicates that net currents flow along the flux tubes. 
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This paper argues that axial uniformity is the result of a rather complex sequence of 
events which occur whenever an electric current / is made to flow along an initially axially 
non-uniform, current-free, axisymmetric magnetic flux tube (a process corresponding to 
injection of magnetic helicity into the flux tube). The sequence of events occurs even when 
/ is modest, i.e., even when the flux tube is only slightly twisted. 
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FIG. 2: (a) Potential (i.e., current-free) coronal loop above solar surface with corresponding sub- 
surface field and source currents. The cross-section of the loop is largest at the top of the loop 
where the magnetic field is weakest; (b) straight cylindrical representation of coronal loop used in 
model. 

The typical arched shape of coronal loops is shown schematically in Fig. 2a; to make 
the analysis tractable we will assume that the loop is straight as sketched in Fig. 2b. 
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However, in order to retain an important aspect of the arched shape, we will allow the 
length of the straight loop to vary to take into account possible variability in the length 
of the arched loop. The Fig. 2b geometry will be characterized by a straight cylindrical 
coordinate system {r, <f), z} where z refers to the direction along the loop axis, is the 
azimuthal direction about the axis, and r is the distance from the axis. The <fi direction is 
called the toroidal direction and the r, z directions are called poloidal. Flux coordinates will 
also be used when appropriate. This poloidal/toroidal nomenclature is formally the same 
as that used for tokamaks, but the configuration should not be confused with a tokamak as 
there are no closed poloidal field lines. The current / will be assumed to be relatively small 
so that the poloidal field magnitude ~ B z is always much larger than in which case the 
flux tube is only slightly twisted. We note that this straight cylindrical approximation of 
flux tube geometry has been used in many previous studies of flux tube equilibria, especially 
force-free equilibria (for example, see Refs. 

BBBBBB) 

. However, the analysis presented 
here differs substantively from these previous studies because our analysis does not begin 
by assuming existence of an equilibrium. Instead, our analysis characterizes the dynamics 
that lead to an equilibrium and shows how the resulting equilibrium is intimately related to 
these dynamics. Furthermore, our analysis takes into account the non-force-free aspects of 
the equilibrium (i.e., finite pressure gradients) and shows that these finite j3 aspects are of 
vital importance to the axial uniformity of the equilibrium. 

Because of the assumed axisymmetry in Fig. 2b, the magnetic field can be expressed as 

B=-^(W>x V0 + /i O /V0) (1) 

where V</> = <p/r and 

rr 

$(r,z,t)= / B z (r',z,t)2Tir'dr' (2) 
Jo 

is the poloidal flux. Axial non-uniformity corresponds to having dip/dz ^ and axial 
bulging corresponds to having ip~ 1 d 2 tp/dz 2 > 0. The current I is similarly given as 

I(r,z,t)= [ J z (r' } z,t)27rr'dr' (3) 
Jo 

and is related to the toroidal field by Ampere's law, i.e., ^1 = litrB^. The current rise time 
t is assumed to be sufficiently slow that Alfven wave propagation effects are unimportant 
i.e., it is assumed that Vat » h where 2h is the length of the flux tube. The current thus 
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flows as in an electric circuit so that there are no retarded time or radiation effects. Taking 
the curl of Eq. ([IJ shows that the current density associated with the magnetic field is 

r 2 S7(b ( 1 \ 1 
J = • -W> + —V/ x V0. (4) 

Controversy exists regarding the properties of I(r, z) external to the current-carrying flux 



;ube. Some authors 1( 



ll|] argue that / must vanish outside the flux tube while others 



12L llH 14, [3] argue that / should be finite outside. If one insists that I vanishes outside 
the flux tube so that there is no net current in a flux tube, the flux tube acts like a coaxial 
cable (i.e., a center conductor sheathed by a coaxial outer conductor carrying equal and 
opposite current). This neutralized current configuration has = external to the flux 
tube and so cannot produce magnetic forces on external currents. Adjacent neutralized flux 
tubes thus cannot exert forces on each other and so will not mutually interact, just like 
adjacent coaxial cables will not mutually interact. Furthermore, sections at different axial 
positions along the length of a neutralized flux tube cannot interact via magnetic forces. 
The lack of interaction between sections at different axial positions along the length of a 
neutralized flux tube means that such a flux tube cannot undergo a kink instability since a 
kink instability involves relaxing to the state of lowest mutual interaction energy between 
loop segments (similarly, a coaxial cable will not undergo a kink instability since there is no 
force between axially separated segments). 

In contrast, having a net current (i.e., being non-neutralized) implies existence of an 
external potential magnetic field ~ J/r outside of the flux tube just like the magnetic 
field external to an ordinary current-carrying wire. An individual flux tube with net current 
can kink since, if its axis bends, there will be forces between sections at different axial 
positions. Two flux tubes, each carrying net current, will experience mutual interaction 
forces due to the B^ of one flux tube interacting with the current of the other flux tube. 
Thus, two adjacent flux tubes each with net current will tend to wrap around each other as 



shown in Fig. 2a of Ref. 



16( since the axis of each flux tube will be affected by the B^ of 



its neighbor. Since we are assuming here that B^ « B z this wrapping will be very slight. 
Examination of the loop structures #1 and #2 denoted by arrows in Fig.l show that these 
two structures do indeed wrap around each other slightly (on the left, loop #1 is to the rear 
of loop #2; on the right, the two loops appear to be in the plane of the photo). 

Based on the observations that kinks do occur in solar structures and that coronal loops 
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do show evidence of wrapping we will assume in this paper that net current does flow in a flux 
tube, i.e., that the current is non- neutralized. This assumption is additionally supported by 
recent work by Feldman Q and by Wheatland H who, after careful analysis of a variety 
of observational evidence, have concluded that net currents do flow. 

The assumption of net current means that we are allowing flux tubes to interact with 
each other via magnetic forces. However, since the effect of this interaction is to alter the 
three dimensional locus of a flux tube axis, and since we are invoking the straight axis ap- 
proximation (i.e. using Fig. 2b to represent Fig. 2a), we are removing from consideration 
the evolution of the three dimensional locus of the flux tube axis. The straight axis approx- 
imation is thus analogous to a kinematics problem where one works in the center of mass 
frame of a body and so removes from consideration external body forces that change the 
location of the center of mass. 

The dynamics of the configuration are governed by the combination of the magnetohy- 
drodynamic (MHD) equation of motion 

p f| =Jx B-VP (5) 

and the induction equation 

°— = V x (U x B) (6) 
where the latter is obtained from the curl of the ideal MHD Ohm's law 

E + UxB=0. (7) 

The flux tube ends are at z = ±h and the flux tube middle is at z = 0. The flux tube 
is assumed to be initially current-free so that initially 1 = 0. The field in the flux tube 
is thus initially potential (i.e., V x B =0 initially) and the source currents generating this 
potential field are external to the flux loop and, as sketched in Fig. 2a, are assumed to be 
below the solar surface. Thus the initially potential coronal loop sketched in Fig. 2a will 
have a magnetic field that is stronger at the footpoints than at the arch top because the 
arch top is further from the source currents. This means that the initial current-free flux 
tube will be bulged at the top because the magnetic field is weaker there. In the straight 
geometry representation given by Fig. 2b, the bulging at the arch top corresponds to having 
the initially potential poloidal field stronger at z = ±/i than at z = and the flux tube 
diameter larger at z = than at z = rkh. 



The sequence of events that occurs when the current is ramped up to a steady-state will 
be shown to consist of the following three stages: 

1. The first stage consists of a twisting of the magnetic field about the z axis in Fig. 
2b together with an associated transient toroidal plasma velocity U^. This stage is 
incompressible and maintains the flux tube poloidal profile, i.e., ip(r,z) is unchanged 
and the flux tube remains bulged. The velocity is proportional to zdl/dt and the 
toroidal acceleration is proportional to zd 2 I/dt 2 . 

2. The second stage involves generation of axial plasma flows. These flows go from z = ±h 
where the flux tube diameter is small to z = where the diameter is large. The flows 
are driven by a ^-directed force which is proportional to —dl 2 /dz. This axial force is 
a nonlinear function of I in contrast to the first stage toroidal acceleration which is a 
linear function of I. 

3. The third stage involves stagnation of the converging flows at z = resulting in plasma 
heating as the flow kinetic energy is converted into heat (thermalized) . There is also 
an accumulation of convected toroidal flux at z = which leads to an enhancement of 
the pinch force at z — 0. The enhanced pinch force squeezes the flux tube diameter at 
z = so that the flux tube approaches axial uniformity, i.e., dip/dz — > 0. Ultimately, 
an axially uniform flux tube loaded with hot plasma results. 



II. LACK OF EQUILIBRIUM FOR ARBITRARILY SPECIFIED MAGNETIC 
FIELDS 

Arbitrarily specified magnetic fields do not, in general, have associated MHD equilibria, 
i.e., in general, no pressure P(r) exists which satisfies 

J x B = VP (8) 

for arbitrarily specified B(r). The essential physics underlying this assertion is that V x 
VP is identically zero whereas V x (J x B) is not necessarily zero, i.e., VP is always a 
conservative force whereas J x B is in general non- conservative. A non-conservative force has 
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an associated torque and since a pressure gradient cannot balance a torque, no equilibrium 
is possible when V x (J x B) is finite. 
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FIG. 3: Poloidal flux surfaces of a magnetic field configuration for which no finite pressure gradient 
MHD equilibrium is possible. 



As a specific example that equilibria do not exist for arbitrarily specified magnetic field 
configurations, consider the simple situation sketched in Fig. 3. We assert that no MHD 
equilibrium is possible for this configuration, i.e., this configuration cannot satisfy Eq. (JSj). 
This assertion is established by calculating radial pressure balance at axial locations z\ and 
Z2 where the respective flux tube radii are a and b with b > a. If satisfied, Eq. (jHJ) imposes 
the requirement B-VP = which means that the pressure must be the same everywhere 
on a field line. 

The poloidal component of the magnetic field is assumed to be straight, uniform, and 
in the z direction at both z\ and Zi . Because the field lines are more densely packed at z\ 
than at Zz, the magnetic field is such that B z (z\) > B z {z<i). In addition, since the poloidal 
magnetic field is uniform at Z\ and Z2, the toroidal current vanishes at both z\ and Z2 [in 
fact, uniformity of poloidal field is more than is needed for the toroidal current to vanish as 
the most general requirement for the toroidal current to vanish is to have V • (r -2 V?/>) = 0]. 

Using = the radial component of Eq. (JHJ) is therefore 

3P/3r = -J.B. = -^L§. (9) 

However, the toroidal component of Eq. (jHJ) gives 

(V/ x V0) x (W x V0) = (10) 

which implies that I = Since the poloidal magnetic field is straight and uniform at 

both Z\ and Z2, the poloidal flux function must have the form ip = ip r 2 /a 2 at Z\ and 
ip = ipor 2 /b 2 at Z2 where ?/>o is the flux on the surface for which P vanishes and a, b are the 
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respective radii of these flux surfaces at z\ and z^- The simplest non-trivial possibility for 
/ = is to assume that / is a linear function of ip so 

HqI = at]). (11) 

Radial pressure balance at z = z\ thus becomes 

dP/dr = = _iM Z_ (12) 

(27ir) 2 fi dr 2ir 2 fi a 4 v ' 

which can be integrated to give the on-axis pressure at Z\ to be 

P(0,zi) = fdrdP/dr = fdrff^ = (13) 
J a Jo 27r 2 /i a 4 4ir 2 a 2 fi 

However, a similar evaluation of the on-axis pressure at z = Z2 gives 

p <°-*> = 43& (14) 

which is smaller than the on-axis pressure at z — Z\ since b > a. Thus, the pressure is not 
uniform along the on-axis field line and so the requirement B • VP = is violated on the 
z-axis. Equation (jHJ) is therefore not satisfied and so equilibrium does not exist. Because 
there is radial pressure balance but not axial pressure balance, we might expect flows to be 
driven from z\ to z-i- The situation is analogous to squeezing the end of a toothpaste tube 
and having toothpaste squirt out the mouth of the tube. 

III. GENERAL CASE USING FLUX COORDINATES 

Let us now return to the problem of a coronal loop which is initially current-free and 
bulging at z = 0, but then has an externally driven current / slowly ramped up to a constant 
value. We define a right-handed orthogonal coordinate system based on the poloidal flux 
coordinates. The unit vectors of this system are 

^ ~ m 

GBpo1 - fvY x V0| (15) 
so that e^p — > f and &Bpoi — ¥ z if the field lines are straight and axial. The form of shows 
that = B • e^, is zero by definition; i.e., the magnetic field can never have a component 
in the direction of V^. 



The current density can be decomposed into the components 



J * = s WxV *-$i (16) 



_ 1 V/ ■ 

We now argue that J^, and Jb po i each have distinctive physics. 

The component flows normal to flux surfaces and provides the torque that causes 
the plasma to rotate toroidally. This current can only be transient and is identified as 
the polarization current Q|. Polarization current can be thought of as being an essentially 
dependent quantity; that is, one first determines the amount of toroidal acceleration using an 
analysis that does not involve the equation of motion, and then one inserts this acceleration 
into the equation of motion to calculate the required polarization current. The reason for 
this inferior status of the toroidal component of the equation of motion is that the toroidal 
symmetry of the system provides a strong constraint on the dynamics. From an MHD point 
of view, the toroidal symmetry means that no toroidal pressure gradient can exist and also 
no toroidal electrostatic electric field can exist. From a particle Hamiltonian point of view, 
this symmetry means that the maximum excursion particles can make from a flux surface 
is no more than a poloidal Larmor radius j^J, a microscopic length. The localization of 
particles to the vicinity of a flux surface means that there cannot be any sustained net 
current density in the direction normal to a flux surface and so is highly constrained. All 
that is allowed is a short-lived transient having zero time-average; i.e., can only be an 
AC current. The slight bobbing back and forth of particles off of a flux surface constitutes 
the polarization current. Thus the plasma acts like a capacitor in the direction normal to the 



flux surfaces, but like a wire in the direction along the flux surfaces. As is well known 21], 
the dielectric constant of the plasma "capacitor" is given by fiopc 2 /B 2 = c 2 /v\ » 1. 
Polarization currents have an associated polarization electric field normal to the flux surface 
resulting from the particles making their small excursions from their nominal flux surface. 

The current / is assumed to be generated by some sub-surface dynamo and so its time- 
dependence is a prescribed quantity. This time dependence is assumed to be such that / 
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increases smoothly from zero to some finite value in a characteristic time r » H/va- This 
smooth increase can be represented by the characteristic time profile 

m = (tanh( t /r) + l) /o (lg) 

so that / = for t « —r and I = I for t » r. It should be noted that dl/dt ~ 
1/ cosh 2 (t/r) which has its maximum at t — and that d 2 I/dt 2 ~ — tanh(t/r) / cosh 2 (t/r) 
which is positive for t slightly before t — and negative for t slightly after t — and then 
otherwise zero. 



Initially untwisted potential flux loop 




Axial current twists flux loop 




FIG. 4: First stage dynamics: Top shows flux tube when 1 = 0. Bottom shows flux tube when / 
is finite, but ip(r, z) is not yet changed from its initial potential value. The evolution from the top 
to the bottom situation involves toroidal rotation by a finite amount, but no poloidal flows. 



A. First stage (ramp-up) 

At the beginning of the first stage / is zero and the flux tube is untwisted as shown in 
the top sketch of Fig. 4. The toroidal component of the induction equation is 

^ = rB por V (^j - rU por V - B^J • V pol (20) 
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and the toroidal component of Ohm's law can be expressed as 

^+U por V^=0. (21) 

We will show that XJ po i results from a weak non-linear force and so is negligible at times 
t ~ because there has not been enough time for a significant XJ po i to develop. In contrast, 
Uff, is proportional to dl/dt and so is large at t ~ when dl/dt is at its maximum. 
However goes to zero at large times when dl/dt — ■> 0. Thus, the time when is finite 
precedes the time when XJ po i is finite. The first stage characterizes the time t ~ when 
Upoi is negligible, if) is unchanged from its initial potential state, is negligible, and U$ is 
transiently finite. 

Letting s denote the distance along the poloidal field from the z = plane and taking 
into account that V po i — 0, Eq. (j20j) may therefore be approximated at t ~ as 

which may be integrated with respect to s to give 

^ _^^> = dl f23l 

The finite toroidal displacement rA<f> = J U^dt is proportional to s and so gives a twisting 
up of the flux tube as shown in the bottom sketch of Fig. 4. In fact this twisting motion is 
such that 

rA(f> s 



(24) 



showing that the plasma twist is just what is required to keep the plasma frozen to the 
twisting magnetic field. 

Equation (J23|) and Eq.(JZJ) together imply the existence of an electric field in the 
direction 

E$ = —U^Bpoi = — — ^ — — ; (25) 
2ixr at 

this is the polarization electric field. The toroidal component of the equation of motion is 

dJT 

~ -J^B pol (26) 

since = 0. Thus, the current normal to a flux surface is 

p dU^ _ p dE^ p_p sd 2 I 

3 pol dt ~ Bl ol dt = Bl ol 



J + B nnl dt BL dt BL27rr dt* • 1 ' 
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Equation (|27|) clearly shows that is indeed the polarization current and that Jf is essen- 
tially a dependent quantity since it is proportional to d 2 I/dt 2 . The polarization current is 
transient and, for positive s, is first negative and then positive (and vice versa for negative 
s). Both Ufj, and the polarization current vanish when / is in steady state. The chain of 
dependence is such that the induction equation first determines which then determines 
Jrf via the equation of motion. For a long, thin flux tube, s ~ z, B po \ ~ B z , and the if) 
direction is approximately the r direction. For t » r, the poloidal current I is in steady 
state and so = for t » r. 



B. Second stage (steady /, development of finite U po i) 



Since I is constant in the second stage, and are both zero. Equation (JTfiJl shows 
that = implies / = in which case surfaces of constant ip(r, z) are also surfaces of 
constant I(r,z). At the beginning of the second stage, XJ po i has not yet developed and so 
ip(r, z) is assumed to be unchanged from its initial potential (vacuum) state, i.e., the poloidal 
profile of the flux surfaces is not yet deformed. Thus, = —r (27r/! ) 1 V ■ (r~ 2 Vip) = at 
the beginning of the second stage. 

We now consider the dynamics. The magnetic force can be decomposed into toroidal and 
poloidal components as follows: 

J X B = J po ;XB po / + J po ^xB tor + J{o r xBp i. (28) 

However, J po ;xB po ; = J^xB po | = and J t0 r = J<t><fi — so that the magnetic force at the 
beginning of the second stage reduces to 

J x B = J pol xB tor = -L (V/ x V<f>) x = -^VJ 2 . (29) 

Since the curl of the magnetic force given in Eq.(|29|) is non-zero, it is impossible for a pressure 
gradient to balance the magnetic force at this stage. The z component of the magnetic force 
is 

which is independent of the sign of /, nonlinear in J, and such as to accelerate plasma from 
regions where the diameter of the current channel is small to regions where the diameter is 
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large. In the case of a flux tube which is bulged in the middle, the force given in Eq. (j30j) 
will accelerate plasma axially from z = ±h towards z — 0. 

The force given in Eq. (J3U|) vanishes at r = since J ~ r 2 for small r. However, axially 
localized radial force balances will quickly develop between the magnetic force and the radial 
pressure gradient as discussed in Section II. The resulting radial force balance will produce 
an axially non-uniform pressure and so there will also be an axial force due to the axial 
pressure gradient. 

Specifically, radial pressure balance means that 

— = -J,B = _^^!. (31) 
dr ^ 8ir 2 r 2 dr 

This radial pressure balance equation is not integrable for arbitrary I(r,z), but to get an 
idea for the general behavior we make the assumption that / ~ r 2 which is integrable. Thus, 
assuming I(r, z) = (r/a(z)) 2 I where Iq is the total current flowing in the flux tube of radius 
a(z), Eq.([3*Tj) can be integrated to give 

so that 

dP _ /i / 2 /I r 2 \ da 



dz ir 2 a 3 \2 a 2 J dz 
is the axial force due to the axial non-uniformity of the pressure. 

Using I(r, z) = I r 2 /a 2 to estimate the axial component of the magnetic force gives 

[J jz " 8vr 2 r 2 dz ~ 2tt 2 a 5 dz [ } 

The total force in the z direction for this case is thus 

dP 



JxB »- Bz 

( r 2 \ da 



2-K 2 a 3 \ a 2 J dz 
This total force is peaked on the axis and has magnitude 

d f\BV 



(35) 



dz 



(36) 



This force will result in axial flows from z = ±h to z = with velocities that are of the order 
of B^r = a, z = ±h) / J^p- Because = 0, the behavior is essentially identical to the 
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situation where tb = and so the flow acceleration mechanism is similar to that discussed 
in Refs. Q-Q which cous.der MHD arc-jets for the situation of purely toro.dal mag uetie 
fields. 




FIG. 5: Evolution of flux tube due to effect of J x B force when / is constant. Solid line shows 
initial constant ip(r,z) flux surface, dashed line shows same flux surface at a later time. Fluid 
elements (shown hatched) are pushed axially towards z = while staying on the same constant tp 
surface. The resulting fluid flows collide at the z = plane where they thermalize their directed 
kinetic energy. Toroidal flux and mass accumulate at the z = layer also. The accumulation of 
toroidal flux at z = increases there and so pinches down the flux tube diameter causing the 
flux tube to become axially uniform. 

C. Third stage (convection of toroidal flux, fluid stagnation, heating, compression) 

The force given by Eq.(|29|) has a finite curl and so cannot be balanced by a pressure 
gradient so long as the original potential profile of i/)(r, z) is maintained. Thus, the only way 
for an equilibrium to develop is for the profile of ip{r, z) to change. This is clearly evident 
from the discussion in Section III(B): Eq. (J35|) shows that axial equilibrium can only occur 
if da/dz — > 0. 

Attainment of equilibrium involves several inter-related hydrodynamic, magnetic, and 
thermodynamic phenomena which are shown schematically in Fig. 5. The solid lines in Fig. 
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5 show a constant ip surface at an early time and the dashed lines show this constant ip 
surface at a later time. Typical fluid elements are shown as hatched parallelograms (cross- 
sections of toroids) and the motion of these elements is seen to consist of both axial and 
radial motion such that each fluid element stays on its own constant ip surface. The non- 
conservative nature of the J x B force is shown in Fig. 5 by the longer length at larger 
\z\ of the arrows representing J x B. The axial motion corresponds to plasma flows which 
ingest plasma at z — ±h, travel towards z — 0, and then converge and stagnate at z = like 
two water jets pointed at each other. The stagnation converts the flow kinetic energy into 
heat and simultaneously increases the plasma density at z = as plasma accumulates there. 
Thus pressure increases at z = 0. Furthermore, toroidal flux embedded in the plasma is 
convected by the axial flows and so there will also be an accumulation of toroidal flux at 
z — 0. This means that the density of the toroidal flux will also increase at z — 0, i.e., 
will increase in the vicinity of z = 0. 

The increase in B^ can be established more rigorously by considering Eq. (|20J) in the 
vicinity of z = and taking into account that (i) — since / is constant, (ii) XJ po i — > 
at z = since the flow stagnates at z = 0, and (iii) V • U po i < near z = since the flows 
are converging at z — 0. Thus, Eq. (j20j) in the vicinity of z = reduces to 

^ ~ -B^V ■ IV (37) 

which shows that B^ must increase since V • XJ po i < (we note that amplification of a 
magnetic field by a converging flow has previously been discussed in Ref.|2q but has not 
otherwise received much attention). In the vicinity of the stagnation layer at z — 0, the 
continuity equation reduces to 

i^ = -V-U poi (38) 
p at 

which can be combined with the induction equation to give 

dB^, ^ B^dp , , 

dt ~ p dt {6J) 

showing that B^ increases in proportion to the increase in mass density at the stagnation 
layer. Since / is constant and 2nrB ( f ) = pol the current channel radius in the vicinity of the 
stagnation layer must decrease as increases to keep rB^ constant. Thus, the bulge of the 
current channel must diminish as sketched in Fig. 5 and, because / = I{ip)-, the bulge of the 
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constant ifi surfaces must also diminish. The result is that the flux tube tends to become 
axially uniform, hot, and dense. 



D. Changes in length 

Making the flux tube axially uniform increases B po i because squeezing the poloidal flux 
surfaces together results in a larger field. Since B po i ~ B z is much larger than B^, it 
would seem that too much energy would have to be invested into squeezing the poloidal flux 
surfaces together. However, if we recall that the loop is really arched and allow the loop 
length to change in such a way that J B po i • dl remains constant where the line integral is 
over the length of the loop, then the loop length 2h will become shorter as B po i increases. 
If f B po i • dl —const., the stored energy in the poloidal field is 

W„nl = J~f B 2 pd dl-ds 



2fio 

1 

2/i 



J B pol ■ dl J . 



B po i ■ ds 

B po i ■ dl 



2^o 

= const. (40) 

It is reasonable to assume that f B po i • dl remains constant, because B po i is produced by 
currents external to the flux tube (e.g., by the subsurface currents sketched in Fig. 2a). 
These source currents may be assumed to stay constant on the time scale during which the 
flux tube undergoes stages 1-3. If one follows the poloidal field along its entire length both 
above and below the solar surface, then it must satisfy Ampere's law 



Hohxt = f B po i ■ dl 



B P oi ■ dl+ / Bpoi ■ dl (41) 

loop J subsfc 

where the contour consists of the loop above the surface and, in addition, the subsurface 
portion; the contour links links the subsurface source current system denoted as I ext . It 
seems reasonable to assume that the subsurface field remains invariant during stages 1-3 
and so f [oop B po i ■ dl must also remain constant. 

Thus, as the poloidal field lines squeeze together to make the flux tube axially uniform, 
the flux tube becomes shorter in such a way as to keep the energy stored in the poloidal field 
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constant. In this manner, no work needs to be done to squeeze the poloidal flux surfaces 
together. One can imagine that the "field line tension" of the poloidal field shortens the 
length of the loop as the poloidal field is made stronger when the poloidal flux surfaces are 
squeezed together. 



E. Ultimate beta 



The z— directed force given by Eq. ljHOj) can be written as 



dz \ /Jo 

The quantity B\j p can be considered as an effective potential energy and so the fact that 
B^/fiQ is large at z = ±/i and small at z = means that there is an effective potential well 
which the plasma falls down as it moves from z = ±h to z — 0. The order of magnitude of the 
resulting flow velocity is given by the reduction in potential energy due to the plasma falling 
down the slopes of the well and so the resulting flow velocity will be U% ~ B^j /x p where B^ is 
evaluated at z = ±h where B 2 is largest. Thus, the flow velocity is of the order of the Alfven 
velocity calculated using the toroidal field (this is much smaller than the Alfven velocity 
calculated using the poloidal field on the assumption that the flux tube is only slightly 
twisted). At the stagnation layer the converging flow velocity is thermalized and so the 
plasma pressure at the stagnation layer will be of the order of P = nmiV^ i +nm e v^ e ~ pU% = 
[Bq/ Ho] . Assuming that B^ « B poh the plasma ft resulting from flow stagnation is 
therefore 

where a is the radius of the flux tube and if) = B po iira 2 . Using the definition a = poI/if> then 

ft = aV/2 (44) 

is the value of ft resulting from flow stagnation. 

The diminishing of the bulge squeezes together the poloidal field so that there will be 
a finite J^, but if the flux tube is squeezed to the point of being axially uniform, then 
vanishes again. Thus starts out by being zero, becomes finite, and then becomes zero 
again if and when the flux tube becomes straight. 
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If an equilibrium is established then J x B = VP which implies B-VP=0 and P = P{ip). 
We define ipo as the flux surface on which P vanishes and tp(r,z) = ip(r,z)/ipo as the 
normalized flux so that 

P(V0 = (l-V0Po (45) 
where Pq is the on-axis pressure. The equilibrium equation J x B = VP can then be written 



in Grad-Shafranov 



27| form as 



r^( 1 -?±) + ^ + a^ = 2 r -l (46) 
dr \r dr J dz 2 a 2 , a 2 , 

where a is the flux tube radius at z = and is defined in terms of the mean axial field at 
z — 0, i.e., = 2/i Po/ (^o/^Oo) 2 • The general solution to Eq. (|47)j) is 

r 2 

iP{r,z) = 2— ^^+ 1X (r,z) (47) 
where x{ r i z ) is an y solution to the homogeneous equation 



r- 



d (ldx \ , 9 2 x 



+ 1^4 + « 2 X = (48) 



9r \r dr J dz 2 

and 7 is a constant to be determined. The condition ip(r, z) = 1 when z = and r = ao 
fixes 7 so that the general solution to the Grad-Shafranov equation is thus 

^( r ' z ) = —TZi + 1 T73 ./„ ( 49 ) 



a 2 ao V a 2 alJ x(ao,0) 

If /3 = a 2 ag/2 then the on/?/ solution to Eq.(|4T)|) satisfying the prescribed boundary condition 
that P vanishes when i[) = i[)q is the particular solution ip = ip^r 2 /a^. However, this solution 
is axially uniform and so we no longer need to specify ao as being the radius at z = 0; it is 
in fact the radius at all z. Equation (|4*9*|) provides the important result that having a finite 
but extremely small will cause the poloidal flux surfaces to differ substantially from the 
force-free situation where is exactly zero and in particular will cause the system to become 
axially uniform when = a 2 a^/2. From a mathematical point of view this is because the 
right hand side of Eq. (j46j) is an inhomogeneous term (source term) of the partial differential 
equation. The source term results in there being a particular solution which would not exist 
if Eq. ()46j) were homogeneous, i.e., if were exactly zero. 

The axial uniformity condition = a 2 a 2 /2 is just Eq. and so we conclude that 

the produced by flow stagnation is precisely the required to force the Grad-Shafranov 
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equation to give an axially uniform solution. We further conclude that, given sufficient time 
and assuming there are no losses, current-carrying flux tubes will always tend to become 



axially uniform and will a 
the well-known property 



ways tend to have the (3 given by Eq. (|44|) . This result is similar to 
28| of tokamaks having (3 po i of order unity because diamagnetism 
exactly balances paramagnetism so that the resulting field is a potential (vacuum) field. The 
roles of poloidal and toroidal directions are interchanged in the coronal loop compared to a 
tokamak and so in the coronal loop it is (3$ which is of order unity. 

The predicted (3 can be compared with actual observed values of (3 in solar coronal loops. 
To make a prediction, a nominal observed flux loop radius a = 1.6 x 10 6 m |2j and a nominal 
measured active region a = 2 x 10~ 8 m _1 are used j^]. These parameters predict a nominal 
(3predicted = ct 2 a 2 /2 = 5 x 10~ 4 . The observed value f3 b serv ^iis calculated using a nominal 
density n = 10 15 m~ 3 , and a nominal temperature 10 6 K |2|. In addition a nominal axial 
magnetic field B z = 1.5 x 10~ 2 T is assumed based on the argument that since the flux 
tube is axially uniform, its axial field must also be axially uniform and so will have the same 
value as the nominal B z = 1.5 x 10~ 2 at the surface of an active region. These parameters 
give ^observed = 2[ionnT / B 2 = 4 x 10~ 4 which is similar to the predicted value. 

This model also has implications regarding the brightening typically observed when the 
axis of a coronal loop starts to writhe and the loop develops a kink instability (sigmoid). 
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32, 



33, 



34j and for a long thin flux tube 



Since kink instability occurs when ah ~ 27r |30, 
a << h, this model predicts that (3 = a 2 a 2 /2 << a 2 h 2 /2 will still be small even if a is 
increased to the point where ah ~ 2tt and kink instability occurs. However, (3 will increase 
as a increases. Since the brightness of a loop is proportional to n 2 for a given temperature, 
this model predicts that the loop should brighten in proportion to the writhing of its axis 
(i.e., in proportion to a as ah approaches unity). 

The model thus provides a heating mechanism (stagnation of MHD-driven flows) which is 
consistent with observed coronal temperatures and densities; however, the prediction is for 
nnT rather than for temperature or density separately; a more detailed analysis would be 
required to isolate the individual dependence of temperature and density on the stagnation 
process. 
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F. Energetic Tail 



As the flows converge, there will be a few particles which have collision mean free paths 
and trajectories such that they bounce back and forth between converging fluid elements. 
Because these particles gain energy on each bounce between the converging flows, these 
particles will gain energy without bound until desynchronized or lost, i.e., they will undergo 
^ acceleration The oa.be, of particles having t ne „a t e mean free path will 
be small, so one will expect a small high energy tail located around z = 0. The concentration 



of 



righ energy particles around z — 0, i.e., at the top of a loop, is in fact what is observed 



IV. SUMMARY AND CONCLUSIONS 



We have shown that the apparently simple problem of driving an electric current along a 
pre-existing potential magnetic flux tube is actually quite complicated and consists of three 
stages. In real situations these stages would overlap and not be as distinct as outlined here. 

The first stage involves a twisting of the magnetic field and an associated ^-dependent 
toroidal rotation (i.e., twisting) of the plasma; this motion is incompressible. The second 
stage involves convergent axial flows driven by the nonlinear, non-conservative force asso- 
ciated with the axial gradient of B^/fiQ. The third stage involves accumulation of both 
mass and toroidal flux at z = and a simultaneous conversion of directed flow energy into 
thermal energy, i.e., stagnation. The concomitant increase in toroidal magnetic field at the 
stagnation layer ultimately leads to an equilibrium and because the flow stagnation gives 
(3 = a 2 a 2 /2 the equilibrium is axially uniform. This sequence of events should be quite 
common and should explain why current-carrying flux tubes are so often observed to be 
filamentary. 

Finally, it should be emphasized that symmetries in both cf) and in z play a critical 
role in the behavior described here. Symmetry in (i.e., toroidal symmetry) prevents the 
existence of any toroidal electrostatic field so that the only allowed toroidal electric field is 
the toroidal electric field associated with changing poloidal flux, i.e., = —(2irr)~ 1 dip/dt. 
Particles are therefore constrained to stay within a poloidal Larmor radius of a flux surface 
so that there can only be AC currents in the direction normal to a flux surface in which case 
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the plasma acts like a capacitor in the direction normal to the flux surfaces. Because toroidal 
acceleration is driven only by the current normal to a flux surface and because no toroidal 
pressure gradient is allowed, the toroidal motion is constrained to be transient, finite, and 
dependent on the temporal behavior of /. Thus when / is constant, poloidal current flows 
along poloidal flux surfaces and there is no toroidal motion. Symmetry about the z = 
plane causes this plane to be a stagnation layer where opposing plasma jets collide resulting 
in accumulation of mass and of frozen-in toroidal magnetic flux and also a thermalization 
of the flow kinetic energy. 

Acknowledgement: Supported by United States Department of Energy Grant DE-FG03- 
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